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Abstract 

The model-independent analysis of the S- and P-wave ttyt phase shifts was carried 
out. This analysis was based on the using of the Roy equations only and all available 
experimental data from the threshold up to dipion mass = 1 GeV. As the results 
S-wave lengths were calculated: Oq = (0.212 ± 0.015) Oq = (—0.043 ± 0.010)m“^. 
The result obtained obviously confirm the standard ChPT version. Moreover, addi¬ 
tional arguments were found in favor of the ratio of the S-wave phase shifts Sq(s) and 
(5 q(s) being independent of energy from the threshold up to m^.,^ = 900 MeV. The 
proportionality coefficient between the phase shifts r] is equal to —4.66 ± 0.05. 
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1 Introduction 


An investigation of the near-threshold parameters of the tttt interaction has acquired a spe¬ 
cial role due to emergence of QCD theories with a broken down chiral symmetry. During 
last years, Chiral Perturbation Theory (ChPT) [H E] and Generalised Chiral Perturbation 
Theory (GChPT) 0 were developed. Both these theories can describe the strong interac¬ 
tions at low energy. The determinative factor in these theories is the existence of vacuum 
condensates violating chiral symmetry. These theories having the same form of the effective 
Lagrangian differ from each other by value of quark condensate and light quark masses. The 
fact determining the choice of the version is that the S-wave tttt scattering lengths aj] and Og 
are very sensitive to the parameters of the model and consequently are the key parameters 
for unambiguous determination of the scenario of chiral symmetry violation. In this way, 
ChPT predicts the value aQ=0.220 and GChPT Og = 0.263^. So, a reliable determination 
of the TTTT lengths enables one to estimate the amount of chiral symmetry violation and to 
choose thereby an adequate version of the theory. 

During some time, despite of large accumulated experimental material on scattering lengths, 
this choice has been difficult to be made. The matter is that the experiment iCe 4 ^ gave 
evidence in favor of GChPT, whereas most vriV — ttitN experiments inclined rather to 
ChPT. 

The aim of our program, begun in [S] and continued in 013, was to choose a true chiral 
version without using additional constraints based on chiral theories. Therefore our calcu¬ 
lation were based on the Roy equations only and all available experimental S- and P-wave 
phase shifts. In our work [H] very large uncertainties of the tttt lengths were obtained that 
prevented making unambiguous choice. In the next paper [3 the additional relation linking 
the S-wave phase shifts was used. This relation was received on the basis of the analysis of 
the S-wave behavior above the threshold only. Theses of chiral theories were not used at all. 
As the result, the accuracy of determination of S-wave lengths Og and Og was considerably 
improved by means of eliminating the correlation between them. The obtained lengths were 
in a good accordance with the standard ChPT version. 

In the present work it will be shown that adding of the new data from the latest KeA E865 jHj 
experiment to the base experimental data set used in 013, makes it possible to improve 
considerably the accuracy of determination of Og and Og and for certain to choose, without 
using additional constraints, the scenario of chiral symmetry violation. 


2 Roy equations 

The using of the general principles of unitarity, analyticity and crossing symmetry is one of 
the seminal approaches to study tttt interaction. For tttt amplitudes, the integral equations 
known as ’’the Roy equations” [9-11] proved to be rather useful on this path. These equations 
determine the real parts of the partial wave amplitudes which satisfy the analyticity and 
crossing symmetry conditions in the —4 < s < 60 range in terms of ttyt amplitude in the 
physical 4 < s < cx)^ region. The Roy equations combined with the unitarity relations 
constitute a system of non-linear singular integral equations. In deriving these equations, 
the dispersion relations with two subtractions at fixed four-momentum transfer t and the 

^The S-wave scattering lengths Oq and Og are given in 

^Here and below, s is the Mandelstam variable, s = 
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crossing symmetry property of the scattering amplitudes were used. In present work the Roy 
equations were solved to get S-wave tttt lengths. All available experimental S- and P-wave 
phase shifts from the threshold up to = 1 GeV were used. And what’s more the new 
high-accuracy data from latest Ke 4 E865 jH] experiment were added to the base experimental 
data set [12-23] used in p. 

For the Sq wave description the phase shifts 5g obtained in the iiN —>■ tittN and iiN —>■ 
ttttA processes [12-17] were adopted. From (T7] the values of the ’’down-flat” set was used 
only. In the region being studied, the ”down-steep” solution coincides with ’’down-flat” one. 
Whereas the ” up-flat” and ” up-steep” versions cannot be described by a smooth curve and 
are strongly differed from the other data used. The results of the Ke 4 HE] experiments 
were used also. 

For the S 2 wave description the phase shifts Sq obtained in the 7r“p —>■ 7r“7r“A’^+ [18-22] 
and TT+p —TT+TT+n (211 processes were adopted. Precise values of the cross sections 
near the threshold were obtained in |221- This permitted to estimate the values Sq in this 
region under the assumption that phase shifts 5 q are known. 

As the result the phase shifts Sq and their uncertainties were calculated by using cross 
sections a^ 7 r(s) and the values of near the threshold from IHEllIl. The resulting values 
are presented in the Table 1 (Appendix A). For the P-wave describing the results obtained 
in the 7r’''7r“ —>■ 7r’''7r“ and —>■ channels [12,14-16] were used. For the case of 

the charged pions, the Roy equations are given by: 

^e//(s) = Af(s) + - /" ^i{x,s)dx + ^i{s), (1) 

^ Ji 

where T/(a;, s) = Jm/g s) + Im fl{x)K 2 i{x, s) + Im fQ[x)Kli[x, s). Explicit expres¬ 

sions for the kernels Kji{x, s) are given in Appendix B. The corrections (pj (s) estimating the 
contributions from the higher waves {I > 2) and from the large mass region were adopted 
from HD. 


(p[](s) = 13 X 10"®(s2 - 16) ± A(p°; A(p° = 5 x 10“®(s^ - 16) 
<Po(s) = 13 X 10“^s(s — 4) ± Apg; Apg = 6 x 10“®s(s — 4) 

The subtraction terms Ag(s) are expressed in terms of the scattering lengths: 


A[](s) = a[] + (2a[] - 5ag); A^(s) = al - (2a[] - 5a^ 
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( 2 ) 

( 3 ) 


We realized the same numerical method to solve the Roy equations as in jHj without using 
iterative procedures. Due to this approach the problem of convergence of the solutions 
was eliminated automatically and the process of calculation of scattering lengths Og and 
Og became absolutely clear. The solution of the Roy equations (HD comprised some steps. 
First, we performed fitting for each phase shift Sj and obtained smooth curves adequately 
describing experimental data. In particular for the S-wave phase shifts expansion © was 
used: 

= A (c;, + ci ,=> + ... + (4) 

A/S 

where q = |a/s — 4 - is c.m. pion momentum and Cl - are free parameters (1=0,2; k = 1-^4). 
We used m=4 because the increase of the number of terms of the series would not improve 
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Figure 1: S-wave tttt phase shifts. The solid curves represent the result of htting in terms of 
expression (4). 


the accuracy of the smoothing. For the So wave, when n = 106 experimental points were 
used, it was obtained: m=4, x^=137.76; m=5, x^=137.28. For the S 2 wave, when n = 28: 
m=4, y;^=36.48; m=5, x^=36.36. Thus, the describing the S-wave phase shifts by means 

of the used polynomial is stable. Experimental values of phase shifts and htting curves iD 
are shown in Fig.l. In the present study, we assume, as in P, that in the energy range 
considered the P-wave is determined by the rho-resonance almost completely. 

On the second stage, the obtained smooth dependencies <5/(s) were used as input for the Roy 
equations m and after integration, the subtraction terms Af (s) were calculated. 

The values of Refo{si) were taken at each experimental point s* where the phase shifts So{si) 
were measured. By solving the Roy equations for each values of Sj, we obtained the values of 
the subtraction terms ho(sj) and their statistical errors cr^/(s) from experimental data on the 
TTTT phase shifts. This errors are determined ultimately by the errors of the phase shifts ho(sj) 
and were calculated by means of the standard rule of propagation of errors. It should be 
emphasized that the expression for the uncertainties c^iis) does not contain the theoretical 
errors A(Pq{s), since they are not, generally speaking, statistical: they change the behavior 
of the function (poi^) simultaneously for all s. Because of this, the theoretical corrections 
(Pg(s) behave as random functions with respect to A(pg(s). Therefore contribution of the 
uncertainties Ap^^s) in the errors of the lengths Oq and Oq should be calculated separately. 
At the conclusion stage we carried out htting of the dependencies Af (sj) using terms (0 
and determined the S-wave mi lengths. Such approach enabled us to study in detail each 
isotopic channel of the Roy equations by evaluating the contribution of each phase shift 
6l in the resulting subtraction terms Ag(s) separately. It is in this way, it was found that 
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Figure 2: Subtraction terms Aq(s) and Ao(s). The straight lines represent the result of htting 
in terms of expression (3). 

the phase shifts 6q, obtained in the ” electronic experiment” [22] , lead to the result which 
contradicts considerably the result obtained by processing the rest of the phase shifts 6q 
data base. Therefore, we did not use the phase shifts from in the present study. This 
problem will be considered below. Now the solutions of the Roy equations obtained for each 
isotopic channel will be given. Hereinafter r is a factor of correlation between Og and Og. In 
the isotopic channel 1=0, it was obtained by htting the subtraction term AQ(sj): 

= 0.207 ±0.009; ag =-0.047 ± 0.005; r = 0.989 (5) 

y;^/NDF=127/106. After taking into consideration the theoretical error A(Pq it was received: 

a° = 0.207 ±0.015; ag =-0.047 ± 0.011; r = 0.980 (6) 

In the isotopic channel 1=2, it was got by htting the subtraction term AQ(sj): 

a° = 0.295 ± 0.042; =-0.022 ± 0.014; r = 0.996 (7) 

y;^/NDF=32.6/25. After taking into consideration the theoretical error A(Pq, it was obtained: 

a[J = 0.295 ± 0.094; ag =-0.022 ± 0.025; r = 0.988 (8) 

The resulting subtraction terms Ag(s) with htting straight lines m are shown in Fig.2. 

The uncertainties of the tttt lengths Og and Ug are dehned both by the statistical errors of 
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values X(^{si), which are expressed through the uncertainties phase shifts by means of the 
standard rule of propagation of errors and by the theoretical uncertainties Aip^^s). 

It may seem that the lengths obtained from the isotopic channel 1=2 are in accordance with 
the GChPT version, but it is not so. The large uncertainties of the obtained lengths do 
not permit to make any choice at all using this channel only. On the other hand, obviously, 
main information about the tttt lengths is contained in the term Aq, because it concentrates 
in itself the overwhelming part of general statistics. We discuss the obtained solutions of 
the Roy equations for each isotopic channel in detail in order to show that these solutions 
((P and (jHI)) are in accordance with each other within the error limits and to demonstrate 
that the subtraction terms Aq(s) are really linear functions of s. For us, it is an additional 
proof that all the calculation steps in the solving the Roy equations and also all preliminary 
work comprising the htting the phase shifts 5/ were carried out correctly. 

The hnal result was obtained on the basis of the use of all available statistics, i.e., by both 
isotopic channels, 1=0 and 1=2: 

= 0.212 ±0.015; =-0.043 ± 0.010 (9) 

with the correlation coefficient r=0.945. Now we can make a preliminary conclusion: the ob¬ 
tained results unambiguously witness in favor of the standard ChPT version and exclude the 
GChPT one. More detailed discussion of the obtained results will be provided in section 4. 
In conclusion of the present section we shall demonstrate the results obtained by using the 
phase shifts from [221 • S' 2 -wave phase shifts are used only from [22] then the following 

is obtained: 

a) 6q - Hoogland [22] only 

/ = 0 : = 0.209 ± 0.015; al = -0.036 ± 0.011; r = 0.978; x^/NDF = 123/106 (10) 

/ = 2 : a° = 0.140 ±0.071; =-0.079 ± 0.016; r = 0.986; /NDF = 32.5/5 (11) 

It is obvious that the results obtained in the different isotopic channels are contradictory 
for the parameter Og. Moreover, a linearity test is not satished - the value of fh® 

channel 1=2 (HH) shows that the subtraction term Ao(s), obtained by using the phase shifts 
Sq from |23j . is not a linear function of s. 

When the both isotopic channels are used then the following is obtained: 

= 0.163 ±0.015; =-0.071 ± 0.009; r = 0.904 (12) 

Thus, when the phase shifts 5g from [22] were used a concordance was absent both between 
the results obtained in the different isotopic channels as well as with the solution (jH)), obtained 
by using the rest of the phase shifts Sq data base, taken from [18-22]. 

The using of the united phase shifts Sq data base from [18-23] does not change the situation 
in principle by force of statistical domination of the phase shifts from [ 22 ] • 

b) Sq - Hoogland |22| ± ah the rest 

/ = 0 : a[] = 0.209 ± 0.015; al = -0.039 ± 0.011; r = 0.979; x^/NDF = 125/106 (13) 
/ = 2 : a° = 0.182 ±0.075; ag =-0.064 ± 0.018; r = 0.992; x^/NDF = 92/32 (14) 

It was got by using both channels: 

a[] = 0.177 ±0.015; ag =-0.064 ± 0.010; r = 0.921 (15) 

It should be noted that adding the phase shifts S^ from [22] in the channel 1=0, which is 
the main source of the information about Og and Ug, leads to the systematic increase of the 
value Og (13). 
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3 Correlation between ^o(s) and ^o(^) 

In the previous section it was shown that the strongly correlated S-wave tttt lengths in the 
result of solution of the Roy equations were obtained. Such correlation, with r = 1, implies 
that the values Oq and Oq are related by a linear dependence. But this fact signihes that 
the phase shifts Sq(s) and 5o('5) ^y force of the near threshold expansion Sq(s) (x a^q must 
be related by a linear dependence too in some energy region near the threshold. We do 
not know only the range of this region . We analyzed the ratio .^(s) = (5o(<s)/(5g(s) for the 
available experimental data to study this problem. No evident dependence on s was found 
in the behavior of ^(s) from the threshold up to s=42, i.e., up to m7r7r=900MeV (Fig.3). 

In such a way, the simplest 0-hypothesis to verify is the hypothesis about proportionality 
phase shifts in some area above the threshold. As the phase shifts 5o(si) and (^o(^j) were 
measured mainly at different energy values, the smoothed curve (Fig.l) representing the 
fitting function (4) was used for calculation of the phase shifts Sq at the points s = Si, where 
the phase shifts were measured. Thus, the ratio of the S-wave phase shifts was calculated 
as = Si). 

The uncertainties were calculated by the standard rule of propagation of errors and £- 
nally they were defined both by errors of phase shifts 5o(si) and was calculated 

by fitting ,^(s) = ? 7 -const, for interval s = 10 -i- 42: rj = —4.66 ± 0.05; x^/NDF=78/82. The 
fitting ,^(s) = ry-const for the interval s = 4-=-42 gives naturally the same value of t], because 
the statistical weights of the points near the threshold are insignificant. In general, large 
uncertainties of the values ^{s) near the threshold (Fig.3) are caused by the fact that the 
phase shifts Sq(s) and 5o('®) have large relative errors in that region. 

So, the proposed 0-hypothesis is conhrmed by means of the statistical proof. And conse¬ 
quently we can conclude that within the present accuracy of the experimental data the ratio 
of S-wave phase shifts does not depend on the energy in the wide enough region. Thus, for 
this energy region, i.e., for s = 10 -i- 42, the relation take place: 

(16) 

where r/ = —4.65 ± 0.05. So, as stated above, from the fact of strong correlation of the tttt 
lengths follows linear dependence of the phase shifts near the threshold . Then we found the 
proportionality between Sq(s) and 5o('®) some region above the threshold: s = 10 -=-42. 
Our only proposal based on these facts is that we deal with the same proportionality, he., 
we believe that the relation (16) is true from the threshold up to s=42. Hence, in force by 
the near threshold expansion (5g (s) oc Og q, the new constraint on scattering lengths follows: 

ao = Val (17) 

So, an opportunity appears using the constraint (17) to eliminate the correlation between 
Og and Og in the process of the subtraction terms fitting. As the result, the solution was 
obtained, which we denote as ”? 7 -solution”: 

a); = 0.211 ±0.005; =-0.0454 ± 0.0010 (18) 

Here it should be emphasized that r^-solution is in accord with the solution (9) obtained 
without using additional constraints. It is a very important point. This signifies that the 
additional constraint (17) relating the tttt lengths does not correct the Roy equations but 
eliminates the correlation only, when subtraction terms are fitting. 
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Figure 3: The ratio of the S-wave phase shifts ^(s) = 5o(s)/hQ(s). The straight line represents 
the constant r] = —4.66. 

In this sense the condition (17) is a new independent constraint on S-wave tttt lengths. We 
stress than the process of obtaining the 77 -solution and the solution (9) is the same right 
up to calculating subtraction terms Aq(s) inclusive. The difference between them consists in 
using the constrain (17) for obtaining the //-solution on the htting step. The solution (9) was 
obtained without using of any additional constraints. The obtained results are presented in 
Fig .6 and Fig. 7. 

4 Discussion and Summary 

Let us analyse obtained results in more detail. We should start by comparison of our re¬ 
sult (9) with the theoretical prediction received in [27], in which ChPT calculations were 
supplemented with the phenomenological representations based on the Roy equations m- 

aK = 0.220 ±0.005; =-0.0444 ± 0.0010 (19) 

These results are in good accord with each other for both parameters Oq within error limits. 
Hence our result (9) certainly witnesses in favor of the standard ChPT version and excludes 
GChPT one, with Oq = 0.263. Thus, the problem of choosing the true ChPT version, in our 
opinion, is solved. 

But it is possible to put a more tough question: whether there is statistically significant 
conformity between the theoretical result (19) and the result of the model-independent 






































Figure 4: The ratio of the S-wave phase shifts .^(s) = 5Q(s)/5g(s) after the filtration, < 0.9. 
The straight line represents the constant r] = —4.66. 



Figure 5: The ratio of the S-wave phase shifts ^{s) = 5q(s)/5q(s) after the hltration, 
< 0.45. The straight line represents the constant rj = —4.66. 
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analysis (9)? It is seen that la contour ellipses do not intersect (Fig.6). Is the hypothesis 
true that these results are statistically consistent one with the other or may be there is a 
signihcant statistical discrepancy of these results? 

Here it is necessary to take into account one feature which distinguish comparison of the 
results in an one-dimensional case and in a plane. In a plane the probability of a random 
variable to get inside a la contour ellipse is equal P=0.39. Certainly it is not enough to 
draw hnal conclusions. Therefore it is more correct to compare 2a contour ellipses. The 
probability to get in such an ellipse is equal P=0.865. The ChPT solution (19) gets in the 
border of the 2a contour ellipse (9) (Fig.6). he., we have to reject the hypothesis that the 
results (9) and (19) are consistent with probability 13.5%. This probability is very large. 
All this taken together forces us to come to a conclusion, that we do not have sufficient base 
to reject a hypothesis about the statistical agreement of the results (9) and (19). Thus, we 
come to the conclusion, that the solutions (9) and (19) are statistically consistent and do 
not contradict each other. 

Let us carry out comparison with other works in which the results of experiment KeA E865 
for calculation of S-wave tttt lengths were used. In the work |H], where the final results 
of this experiment were presented, it was received without using of the additional relations 
linking Oq and a^: 

a° = 0.203 ± 0.033 ± 0.004,^^,*; = -0.055 ± 0.023 ± 0.003,y,t (20) 

he., we have full conformity with our result (9) within the limits of errors (Fig. 7). In the 
work [2Sj the position and the borders of the area in the plane (oq, Oq) in which S-wave 
lengths are consistent with the Roy equations solution and the available experimental data 
on TTTT phase shifts above 0.8 GeV were specified. It was received for the central curve of this 
area: 

al = -0.0849 + 0.232a° - 0.0865(0°)^ [±0.0088] (21) 

The value given in brackets defines the width of the band. In Fig.6 and Fig.7 this band is 
designated as UB (universal band). 

In the work [22] the calculations done in were repeated with some changes and practically 
the same parameters describing UB were received. Further, using the obtained parametriza- 
tion and the experimental data including the data [Sj, the authors received: 

a[! = 0.228 ± 0.013; =-0.0380 ± 0.0044 (22) 

with the factor of correlation r=0.799. In the work [Sj using UB [22] as the additional 
constraint close results were received. 

The results obtained in [H| and [221 are given in Fig.7. The solution (22) gets in our 2a 
contour ellipse as well as our solution (9) gets in 2a contour ellipse of the solution (22). 
Thus, it is possible to state that the results (9) and (22) do not contradict one another. 

Let us consider the problem of stability of the received solution (9) concerning the procedure 
of experimental data selection, i.e., tttt phase shifts, which in our method of the solution 
of the Roy equations are utilized as input. Stability of the solutions versus variations of 
the initial data is an important indicator of reliability of the method of the solution and 
consistency of the initial data. We have shown above that the use of the data from PI leads 
to contradictious results. Further, the results of an expanded analysis are presented. 
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Figure 6 : (Color online) The S-wave ttyt lengths. The olive ellipse with the centre as a square 
indicates the solution (9); the solid line - Icr ellipse, the dashed one - 2a. The red ellipse with 
the centre as a circle represents 77 -solution (18). The blue ellipse with the centre as a cross 
shows the ChPT result j2Zj; the solid line - la ellipse, the dashed one - 2a. The straight 

lines marked UB indicates the area allowed for S-wave lengths The strip marked CLG 

is the range corresponding the chiral constraint m- 

1 ) Change of the data sets used. 

a) The solution of the Roy equations without the phase shifts Sq which calculated on the 
basis of the cross sections received in [21] (Tabl. 1): 

= 0.213 ±0.015; ag =-0.044 ± 0.011; r = 0.961 (23) 

b) Calculation of S-wave lengths without the data from the work jHj: 

a)J = 0.208 ± 0.015; ag =-0.045 ± 0.010; r = 0.944 (24) 

2 ) Change of the degree of the htting S-wave phase shifts polynomial (4): 

For m=5 in the formula (4) we receive: 

= 0.211 ±0.013; =-0.044 ± 0.009; r = 0.918 (25) 

Comparison of the results (23-25) with the above solution (9) shows, that the criterion of 
the stability for the given solution is satisfied. 

Let us proceed to the discussion of observable proportionality of S-wave phase shifts. It may 
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Figure 7: (Color online) The S-wave ttyt lengths. The olive ellipse with the centre as a square 
indicates the solution (9); the solid line - Icr ellipse, the dashed one - 2a. The black ellipse 
with the centre as a rhomb represents the result [211 • The red dotted ellipse with the centre 
as a circle represents the result jH]. The blue ellipse with the centre as a cross shows the 
ChPT result j2Zj. 

seem, that the values ^(sj) in Fig.3 have a wide scatter and, therefore, can be described not 
only by a constant, but also by some class of smooth functions of s. But these doubts are 
based on visual illusion. The point is that the values ^(sj) with the large errors a^ form ”a 
cloud” which masks true dependence. These points have small statistical weights and do not 
give contribution to r] value. We have carried out a hltration leaving only the points with 
errors less than given, i.e., with a^ < a^ where ak lay in the range 1 0.3. The remaining 

after the hltration values of were htted by a constant. The result: all received lay 
in the range (—4.66 -i—4.68) and have errors cx^ = 0.05 and goodness-of-ht test is satished: 

< iVfc, where Nk - number of points ^(s*) after the hltration with parameter ak- he., all 
sets of points after the hltration are well described by a constant. The results are given in 
Fig. 3-5. Fig. 3 - the values ^(sj) without the hltration, Fig.4 - ak = 0.9, Fig.5 - ak = 0.45. 
One can easily see that after rejection of the points with large errors, the remaining points 
more and more concentrate near the straight line. Thus, if the 0-hypothesis is that values 
.^(sj) are the constant within the considered region, this hypothesis are proved both statis¬ 
tically and visually. 

Additionally, a linear function was utilized for htting .^(sj) also to hnd any dependence of ^ 
on s, if it exists nevertheless. The linear function may be represented as f{s = rj' + b{s — 4), 
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where rj' and b are free parameters. It was obtained as the resnlt of htting for s = 10 -h 42: 

V =-4.51 ±0.15; 6 =-0.006 ±0.006 (26) 

X^/NDF=77/81. As | 6 | < cJb, there is no reason to believe, that the hypothesis abont the 
linear dependence ^(s) is confirmed. 

Let ns consider the problem of inflnence of the rj errors on the ttyt lengths errors obtained by 
the solving of the Roy eqnations (//-solntion (18)). It may seem that the sizes of Oq errors are 
small becanse the sizes are small. Bnt it is not so. The basic contribntion to errors of the 
valnes and Og, in this case, is bronght by theoretical nncertainties in the Roy eqnations 
A(Pq{s). Role of the additional constraint (17) is only to eliminate the correlation between 
Oq and Og in the process of htting of snbtraction terms Aq(s). To show it we increase fonr 
times, i.e., we nsed cr^ = 0.2. In resnlt the following 77 -solntion is obtained: 

aK = 0.211 ±0.0052; =-0.0454 ± 0.0016 (27) 

The solntion (27) shows a weak dependence on 

In the works 123121 it was shown that the width of the allowed area in the (og, Og) plane can 
be rednced considerably by nsing the additional chiral constraint imposed on S-wave lengths. 
This constraint links the combination (2ag — hog) with the scalar pion radins < r^s >• In the 
resnlt of ntilizing of this constraint it was received: 

Aa^ = 0.236Aa° - 0.61(Aa°)2 - 9.9(Aa°)^ [±0.0008] (28) 

where Aoq = Og — 0.22; Aog = Og ± 0.0444. In Fig .6 this narrow strip is denoted as CLG. 
From this hgure it follows that r^-solntion (18) lay practically in the border of the CLG-band 
and half of la contonr ellipse overlaps this band. Also r^-solntion lays practically in the bor¬ 
der of 2a contonr ellipse of the GhPT-solntion (19). Owing to all aforementioned, one may 
conclnde that //-solntion (18) received nnder the additional condition (17) is consistent both 
with the chiral GLG constraint (28) within the la level and with the GhPT-solntion (19) 
within the 2a level . 

It is natnral to compare the value of // = —4.66 ± 0.05 received in the present work with 
f]chPT = flo/flo) which follows from the chiral theory. Only it is necessary to take into account 
that the calculation of the ratio of S-wave /ftt lengths should be carried out in view of their 
correlation. From values Og, Og and 2aQ — hog obtained within GhPT framework [2111201 , one 
may estimate the factor of correlation between Og and Og. It is equal 0.8. In view of it for 
the S-wave tttt lengths ratio it was received rjchPT = —4.95 ± 0 . 21 . The difference from the 
T] value received by us is slightly more than one a. 

Summarizing the main results of the present study, one may say that the solutions received 
by us and other authors [3 1201 grouped near the GhPT-solntion m and are 

consistent both with each other and with this GhPT solution. Thus, we believe that the 
problem of choosing of the scenario of chiral symmetry violation is solved. The available 
mismatch on the la level both among the considered Roy equations solutions and among 
these solutions and the theoretical prediction 1271 may be caused by the fact that we used 
non identical sets of the experimental data and different methods of the Roy equations solu¬ 
tion. Therefore, it seems that prior to search for the physical reasons of such divergence, it 
is necessary to come to an agreement about using of uniform experimental data base. Also 
it is desirable to organize the procedure of the Roy equations solution in such a way that 
enables to check both individual solutions in every isotopic channel and monitor influence 
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of various errors (statistical, systematic, theoretical, errors from additional constraints) on 
the resulting errors of the S-wave tttt lengths. May be that such unihcation of the initial 
data and more detailed control of the course of the solution will allow to reduce the existing 
discrepancy. 

This work was supported in part by the Russian Foundation for Basic Research (project no. 
00-02-17852). 

Appendix A 


Table 1. 


s 

E, MeV 

(5g, deg 

deg 

4.15 

284.3 

-1.11 

0.67 

4.45 

294.3 

-1.75 

0.96 

4.75 

304.2 

-1.05 

1.11 

5.05 

313.7 

-0.71 

1.26 

5.35 

322.8 

0.93 

2.30 

5.65 

331.8 

-5.20 

3.27 
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